Plasma-based generation of a single few-cycle, high-energy and 
ultrahigh intensity laser pulse 
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A laser-boosted relativistic solid-density paraboloidal foil is known to efficiently reflect and focus a 
counterpropagating laser pulse. We show that in the case of an ultrarelativistic counterpropagating 
pulse, a high-energy and ultrahigh intensity reflected pulse can be more effectively generated by 
a relatively slow and heavy foil than by a fast and light one. This counterintuitive feature is 
explained with the larger reflectivity of a heavy foil, which compensates for the lower relativistic 
Doppler factor. Moreover, since the counterpropagating pulse is ultrarelativistic, the foil is abruptly 
dispersed and only the first few cycles of the counterpropagating pulse are reflected. Our three- 
dimensional particle-in-cell simulations show that a single sub-5-femtosecond, multi-petawatt laser 
pulse with several joule of energy and with peak intensity exceeding 10 24 W cm" 2 can be generated 
employing laser pulses with peak intensity of the order of 10 22 W cm" 2 . In addition, the carrier 
envelope phase of the generated pulse can be tuned provided that the carrier envelope phase of the 
initial counterpropagating pulse is controlled. 

PACS numbers: 52.59.Ye, 52.38.-r, 52.65. Rr, 42.65.Re 
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A wide range of novel studies in nonlinear optics as 
well as the major new regimes of extreme field physics 
require laser pulses with three key features: they should 
be few-cycle, high-energy and of ultrahigh intensity. Al- 
ready in the realm of non-relativistic atomic and molec- 
ular physics, the electron dynamics in the presence of 
a few-cycle laser pulse differs qualitatively from that 
in a quasi-monochromatic pulse, and it becomes sen- 
sitive to the precise shape of the driving electromag- 
netic field In turn, this opens the possibility of 
controlling atomic quantum processes by manipulating 
the laser pulse shape. In order to achieve the same 
goal also in the ultrarelativistic regime and in the realm 
of nonlinear quantum electrodynamics, few-cycle laser 
pulses are required with peak intensities largely exceed- 
ing 10 20 W/cm 2 0]. At such high intensities, the spin 
dynamics of an electron in a few-cycle laser pulse is pre- 
dicted to be altered by the electron self- field Q, and 
the nonlinear Compton emission spectrum is expected to 
show pronounced pulse-shape effects Moreover, the 
electron-positron photo-production rate in a few-cycle 
laser pulse may exceed by orders of magnitude that in 
a monochromatic laser pulse 0, and the energy spec- 
trum of the emitted pairs becomes sensitive to the laser 
pulse shape 0- In a more indirect way, few-cycle, high- 
energy and relativistically intense laser pulses interact- 
ing with solid-density plasma surfaces provide a reliable 
and efficient route to the generation of single attosecond 
pulses with unprecedented intensities [l], Q, paving the 
way to the investigation of, e.g., non-pcrturbative and 
relativistic inner-shell processes, and to the development 
of nonlinear XUV attosecond metrology [|| . 

Present Thsapphirc laser technology provides laser 
pulses with peak power up to the PW scale and 30 fs 



duration (loj . However, the limited bandwidth render 
the generation of few-cycle pulses very difficult. Opti- 
cal parametric chirped-pulse amplification circumvents 
this limitation and sub-10-fs duration laser pulses with 
comparatively low energy can be generated [ill ]. The 



generation of few-cycle laser pulses with multi-joule en- 
ergy remains especially challenging. Relativistic flying 
mirrors interacting with counterpropagating laser pulses 
have been proposed in a variety of different configurations 
as a promising frequency multiplication method to pro- 
duce X-ray and gamma-ray beams IT3. HI , to generate 
ultrawide- bandwidth chirped pulses [14 1 . nearly single- 
cycle low-energy pulses |l5j and to attain ultrahigh in- 
tensity, yet low-energy laser pulses by focusing the par- 
tially reflected pulse to its diffraction limit However, 
in these studies the emphasis was on frequency multipli- 
cation and spectral broadening, and comparatively light 
plasma mirrors interacting with relatively low-intensity 
counterpropagating pulses have therefore been consid- 
ered. In addition, in the proposed schemes the mirror 
is only slightly perturbed by the low-intensity counter- 
propagating laser pulse 0, 03 ■ Hence, the energy trans- 
fer from the relativistic mirror to the counterpropagat- 
ing pulse is modest and the energy of the reflected pulse 
remains relatively small. We mention that, employing a 
method based on the ionization-induced self-compression 
effect, the generation of a 6 fs, 0.3-joule pulse has been re- 
cently demonstrated numerically [17(. To our knowledge, 
none of the proposed variants have hitherto achieved 
the generation of few-cycle pulses with multi-PW power, 
multi-joule energy and intensities beyond 10 24 W cm -2 . 

In this Letter, we put forward the concept of a laser- 
boosted solid-density paraboloidal relativistic mirror, in- 
teracting with a superintense counterpropagating laser 
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FIG. 1. Scheme of the proposed setup. 



pulse, to generate a sub-three-cycle pulse with multi- 
joule-energy and peak intensity exceeding 10 24 W cm -2 . 
It is found that in order to maximize the peak intensity 
and the energy of the reflected pulse it is convenient to 
employ a heavy and therefore relatively slow "mirror" . 
This counterintuitive result is explained with the larger 
reflectivity of a heavier "mirror" , which compensates for 
its smaller velocity. In addition, the short time dura- 
tion of the reflected pulse is achieved by employing a su- 
perintense, high-contrast incident pulse, which abruptly 
disperses the plasma mirror after only the first few cy- 
cles. Two-dimensional (2D) and three-dimensional (3D) 
particle-in-cell (PIC) simulations indicate the feasibility 
of the presented setup already at incident intensities of 
the order of a few times 10 22 W cm~ 2 . 

Figure [T] shows the scheme of the proposed setup, 
where a "driver" pulse accelerates a "mirror" along the 
positive x-direction and a "reflected" pulse is generated 
in the collision of the mirror with a counterpropagating 
"source" pulse. Hereafter, the subscript s (d) and the up- 
per (lower) sign refer to the source (driver) pulse counter- 
propagating (copropagating) with respect to the mirror. 
If a laser pulse of frequency ui s /d and intensity I a u is 
normally incident on a flat foil with reflectivity 7Z s /d and 
moving along the x-axis with velocity /3 = v x /c > 0, the 
frequency u> r and the intensity I r of the reflected pulse 
are 
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respectively, with D = + — P) > 1 being the 

Doppler factor. Hereafter, we assume that the source and 
the driver pulse have the same frequency uj s = u c i = lo. 
For a counterpropagating pulse the reflectivity of the foil 
decreases at increasing foil velocities and it also depends 
on the pulse intensity. Our aim here is to determine 
the conditions for maximizing both the intensity and the 
energy of the reflected pulse. In order to develop an ana- 
lytical model, we describe a foil with electron density n e 
and with thickness £ by a 5-like density profile n(x) = 
n e £6(x) ljl, 2(J. The reflectivity of a (5-like foil moving 



with velocity f3 and for an impinging pulse of relativis- 
tic intensity, i.e. with a 2 , s , d = 4ire 2 I s /d/m 2 uj 2 c 3 ^> I, 
reads El 



r 2 
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and Cs/d = Co/-D ±1 - Here, (o is the normalized surface 
density (q = lo 2 £/2luc, with to 2 = / lire 2 n e /m e being the 
square of the plasma frequency. Since it is understood 
that both the source and the driver fields are ultrarela- 
tivistic (a 2 , s , d S> 1), Eq. ((2]) can be approximated as [20j [: 
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Eq. ((3]) shows that in order to efficiently reflect an ul- 
trarelativistic pulse it is convenient to have a relatively 
"heavy" foil with ( s/d > a ,s/d- 

In our model the foil is initially at rest and it is accel- 
erated along the positive x-direction by the driver pulse. 
In order to determine the value of the Doppler factor af- 
ter the acceleration phase Dq as a function of the driver 
pulse parameters, we assume that (d > ao,rf an d thus 
TZd ~ 1 (sec Eq. (|4])). The velocity of a foil accelerated 
by the radiation pressure [2l| of the driver pulse can be 
calculated analytically by employing the "light sail" (LS) 
equation for a perfectly reflecting mirror [iol [22| : 



dp _ 2-irZme 11-/3 a 2 d [w(t)} 
~dt ~ Am p ^l + P <o 



(5) 



where 7 = 1/ y^l — /3 2 and ad[w(t)] is the field amplitude 
of the driver pulse as a function of the foil phase w(t) = 
t — x{t). In Eq. (|5|), Z and A are the ion atomic num- 
ber and atomic weight respectively, and normalized units 
are used: time is in units of the laser period T = 2tt/ui, 
space in units of the laser wavelength A, velocities are in 
units of c and fields in units of m e ujc/\e\. Eq. ([S]) can be 
solved analytically (see e.g. Ref. |20[) and it is found that 
Do = 1 + £d/(o, where Ed = 2irZm e f a 2 d (u) du / Am p is 
the 'effective' energy of the driver pulse. For a sharp-rise 
source pulse, the peak intensity of the reflected pulse is 
reached close to the peak intensity of the source pulse. 
Since during the source pulse-foil interaction the foil un- 
dergoes a strong recoil due to the radiation pressure of 
the source pulse, the Doppler factor D of the foil at the 
peak of the source pulse is smaller than Dq, consequently 
reducing the maximal intensity of the reflected pulse. On 
this respect and also to avoid a premature foil dispersion, 
it is convenient to employ a sharp-rising, high-contrast 
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source pulse, as those generated with the plasma mirror 
technique 23j, |2J| . By proceeding as for the calculation 
of D , we obtain 
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where £ s = 2itZm e j a 2 (u) du/Am p is that part of the 
source pulse energy that is before the source pulse peak 
intensity. Since we seek 1Z S 1, we require that Cs > 
do si which provides the constraint Co > Co.m with 



Hence, lower peak intensity values compared to the model 
predictions are attained. In the considered range of in- 
tensities (I s < 10 23 Wcm -2 ) the foil is nonrelativistic 
during the acceleration stage. In the nonrelativistic case 
(6) the instability grows more quickly for increasing a 2 , d /Co 
(in the ultrarelativistic case more slowly for increasing 
a o d/Co, see Ref. (25)]). Our 2D simulations show that via 
increasing by about 15% the value of Co compared to the 
'optimal' value Co,m, the effect of the RT instability is 
sufficiently reduced while retaining a noticeable intensity 
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where the conditions £d, ao,s 3> £ s have been assumed. 
In order to maximize the reflected intensity at 1Z S w 1, 
we have to maximize the Doppler factor D as a func- 
tion of Co for given £4 and £ s . At £d,an tS 3> £ s the 
value Co at which D is maximal is always less than Co,m- 
Since the Doppler factor D monotonically decreases for 
Co > Co ( see Eq. (HJ)), its maximal value D m compat- 
ible with the condition Co > Co.m is at Co,m, and it is 
D m = Co,m/io.s- The fact that there exists an optimal 
value of the normalized surface density has a simple phys- 
ical interpretation: for fixed driver and source pulses, if 
Co is too large, the foil slows down and the Doppler factor 
is small. If Co becomes too small, the velocity of the foil 
increases and the reflectivity rapidly decreases because 
C s tends to vanish (see Eq. (J2J). Finally, we note that 
since the foil is at rest before the driver pulse impinges 
on it, the condition IZd « 1 can be fulfilled at Co > a o,d- 
Hence, the condition Co,m > ao.d has to be secured, which 
provides a 0jS > a od/( a o.d + £<i)- in this region of pa- 
rameters, the maximal reflected intensity for a flat foil 
l r ^ m oc D m a^ s is an increasing function of I s cx a 2 , s . 
In a 3D geometry, a paraboloidal foil can focus a laser 
pulse to its diffraction limit. Since the laser wavelength 
is Doppler reduced in the rest frame of the foil, the sur- 
face of the diffraction limited focal spot is reduced by 
a factor D 2 . Hence, the maximal focused intensity be- 
comes I r ,m oc D^clq s , which is an increasing function of 
I s oc a 2 , s for a 0iS > £<i/2. In both cases, the driver and 
the source pulse should be superintense, short and with a 
sharp intensity rise in order to increase the peak intensity 
and the energy of the reflected pulse as well as to reduce 
the effect of instabilities [25| . 

In order to validate the model, test the scalability of 
the scheme and check the convergence of numerical re- 
sults, we started with one-dimensional (ID) and 2D PIC 
simulations as they are computationally less expensive. 
Our ID simulation results are in good agreement with the 
above analytical predictions for the peak intensity of the 
reflected pulse. In 2D simulations the onset of Rayleigh- 
Taylor-like (RT) [2l| instabilities lead to the formation 
of high-density clumps separated by lower density re- 
gions 25J, rendering the foil 'porous' to the source pulse. 



increase. 

To confirm the reliability of our 2D simulation results, 
we performed fully 3D simulations with analogous pa- 
rameters. In our simulations the driver laser pulse is cir- 
cularly polarized (CP), its field amplitude has a 'trape- 
zoidal' shape in time with one-cycle sin 2 -function rise 
and fall, and with a five-cycle constant plateau. Its 
transverse field profile at the focal plane is Gaussian 
with spot radius ad = 4.25A, wavelength A = 0.8 /im 
and peak intensity I d = 7.3 x 10 22 W cm~ 2 . The foil 
consists of fully ionized carbon with electron density 
n e = 400 n c , with n c = m e uj 2 /4we 2 being the critical 
density, which corresponds to the density of graphite 
p ~ 2.3 gem -3 . The foil is initially localized at x = 3A 
and it is shaped transversely with a thickness distribution 
£ = max[^i,£ exp(-r 2 /2cr 2 )] with l x = 0.02A, £ = 0.2A, 

(Tf = 2.43A and r = \/y 2 + z 2 . We point out that the 
structure and properties of such nanometer-scale thick- 
ness carbon foil can be engineered with high precision 
with present-day technology [2f||22|- For <Jd > V%<Jf the 
acceleration factor [28( a 2 (f)/Co( r ) ^ s l ar g er m the outer 
part of the foil, which therefore takes a focusing pro- 
file for the counterpropagating source pulse. The effec- 
tiveness and scalability of this approach were confirmed 
with 2D simulations for various laser and foil parameters 
and down to ad = 2.28A and aj = 1.5A. The size of 
the computational box is 15A(x) x 18A(y) x 18A(2) and 
the corresponding grid is 3000(x) x 360(y) x 360(2). We 
use 729 particles-per-cell for each species (electrons and 
ions). The source pulse intensity as well as its temporal 
and transverse profile are the same as the driver pulse 
except that the plateau is three-cycle and its width pa- 
rameter is chosen to match the foil spot a s = aj. The 
driver laser-pulse front reaches the edge of the foil at 
t = 0, while the source pulse reaches the foil at nearly 
the end of the acceleration phase, which lasts approxi- 
mately 8 T. An accurate synchronization between two 
laser pulses can be achieved, e.g., by generating the two 
pulses from the same seed pulse before the amplification 
stage. From the above-mentioned parameters we have 
ao, s = a>o,d — 130, £ s ~ 10.8 and £d — 155.4 such that 
the conditions £d,ao,s ^> £ s assumed for Eq. ([7|) are ful- 
filled. From Eq. we have Co,™ ^ 221.3 while the 
normalized surface density is Co — 251.3, i.e. approxi- 
mately 14% larger than Co,m- Hence, from Eq. ^ we 
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FIG. 2. 3D PIC simulation results showing the reflected pulse 
at t = 14.71 T in the case of a LP source pulse with polariza- 
tion along the j/-axis. The reflected pulse propagates along 
the positive s-axis, see the main text for details on param- 
eters. Panel (a): 3D distribution of the square root of the 
energy density £ r = W (E}. + Bp)/2 in normalized units. The 
energy density distribution is not plotted in the region y < 
and 2 > to show its values inside the reflected pulse. Panel 
(b): j/-component of the electric E r , y (black) and superim- 
posed the z-component of the magnetic B r , z field (red) along 
the central axis. Panel (c): modulus of the Fourier transform 
of the y-component of the electric field \E r , y (k x )\. Panel (d): 
the power P r in units of PW contained in a square with sides 
of length 2A centered on the axis. 



obtain D ~ 1.51 and assuming an effective reflecting sur- 
face S ps 7ro-g ~ 4-7rA 2 the estimated frequency and peak 
intensity of the reflected pulse are: to r = D 2 oj ~ 2.28w 
and I r ps D 6 I S S/X 2 ~ 10 25 W cm" 2 , respectively. 

Figure [2] shows the energy density distribution, tem- 
poral profile, spectrum and power of the reflected pulse 
at t = 14.71 T obtained from 3D PIC simulations. The 
initial source pulse is linearly polarized (LP) along the 
y-axis. Fig. [21a) displays the square root of the energy 
density distribution £ r = (E 2 + B 2 )/2 in normalized 
units. The focusing of the laser pulse close to its diffrac- 
tion limit is evident, the peak intensity of the reflected 
pulse is I r ~ 3.3 x 10 24 W cm" 2 (I r ~ 2.4 x 10 24 W cm" 2 
in the case of a CP source pulse). Fig. [2fb) shows the 
y-componcnt of the electric field E r , y (black) and the 
z-component of the magnetic field B TyZ (red) along the 
central axis. The temporal duration of the main peak is 
approximately 2.7 fs and is followed by a tail of compa- 
rable duration. It is shown that both the main pulse and 
the pulse tail contain only two cycles. The modulus of 
the Fourier transform of the y-componcnt of the electric 
field along the axis \E rjV (k x )\, where k x denotes the wave 
number and k = 2ir/\, is reported in Fig. Etc) showing 
that the reflected pulse is chirped with a remarkable in- 
crease of the spectral bandwidth that spans from the near 



infrared to the far ultraviolet with the main peak at the 
boundary between the visible and the ultraviolet region 
(~ 370 nm) . Both the temporal profile and the spectrum 
are similar for the CP source pulse. Fig. EJd) shows the 
power P r contained in a square with sides of length 2A 
centered on the axis. The sharp peak corresponding to 
the reflected pulse is clearly visible, the peak power be- 
ing 5 PW (5.5 PW for the CP source pulse). Finally, the 
reflected pulse energy is 4.1 J for the LP source pulse and 
7.4 J for the CP source pulse. Comparing the 3D PIC 
simulation results with the model predictions it follows 
that the prediction of the reflected frequency is in good 
agreement with the main peak of the reflected pulse (see 
Fig- EJc)). On the other hand, the intensity enhance- 
ment factor of the reflected pulse is three times smaller 
for the LP source pulse (four times for the CP source 
pulse) than predicted by the analytical model. This may 
be explained noticing that, due to transverse instabilities, 
the foil tends to be torn into higher-density clumps sep- 
arated by lower density regions. Part of the source pulse 
is either transmitted through the clumps or absorbed by 
the foil therefore reducing the intensity and the energy 
of the reflected pulse. 

We stress that in our simulations the physical features 
of the reflected pulse such as the peak intensity, energy 
and number of cycles are not affected by the duration of 
the source pulse, e.g., no appreciable changes are found 
doubling the duration of the source pulse plateau. In 
addition, the reflected pulse 'inherits' the carrier enve- 
lope phase (CEP) of the source pulse. Since the reflected 
pulse contains only the first cycles of the source pulse, 
its physical properties significantly depend on the CEP 
of the source pulse. This also implies that the CEP of 
the reflected pulse can be tuned, provided that the CEP 
of the source pulse is controlled. 

Following Refs. [29|, [3(j, the role of radiation reaction 
(RR) effects was checked comparing the results of simu- 
lations with and without the Landau-Lifshitz force [3l| . 
Simulations show that RR effects do not lead to signifi- 
cant differences in the reflected pulse both for a CP and 
for a LP source pulse. This may be explained noticing 
that when the reflected pulse is generated, the foil den- 
sity is still high such that the skin-depth S ~ c/u p is 
small and the fields inside the foil are much smaller than 
the vacuum fields [2^]. Simulations show that the sole 
signature of RR effects in the reflected pulse is a modest 
intensity increase on the reflected pulse tail,_probably due 
to the 'bunching' property of the RR force [32|, which in- 
creases the electron density and therefore the reflectivity 
of the foil remnants. For the parameters of the simula- 
tions reported here, we ensured that the probability of 
electron-positron pair production is negligibly small. 

In conclusion, a realizable route to the generation of 
sub-5-femtosecond, few-cycle laser pulses with several- 
joule-energy, multi-PW peak power and peak intensities 
exceeding 10 24 W cm -2 has been proposed. Laser pulses 
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with such extreme properties may open new horizons in 
attosecond physics and may allow to probe unexplored 
regimes of nonlinear quantum electrodynamics. Our mul- 
tidimensional simulations show that the proposed scheme 
is testable with next-generation lasers. 

We acknowledge useful discussions with 
B. M. Hegelich, N. Kumar, A. Macchi and G. Sarri. 
PIC simulations were performed using the computing 
resources granted by the VSR of the Research Center 
Julich under the project HRO01. 



[1] 

[2] 

[3 

[5 
[6 
V. 
[8 
[9 
[10 

I" 
[12 
[13 



matteo.tamburini@mpi-hd.mpg.de 
F. Krausz and 



M. 



|Rev. Mod. Phys. 81, 163 (2009) 



Ivanov, 



A. Di Piazza, C. Miiller, K. Z. Hatsagortsyan, and C. H. 
Keitel, Rev. Mod. Phys. in press. ArXiv e-prints (2011), 
|arXiv:1111.3886 [hep-ph]j 

s! Meuren and A. Di Piazza, 

|Phys. Rev. Lett. 107, 260401 (2011)' 



Di 



Piazza, 



V. 



Florescu, 



and A. Hosaka, 



R Mackenroth and A~ 
|Phys. Rev. A 83, 032106 (2011)| 
M. Boca, V. Dinu, and 
|Phys. Rev. A 86, 013414 (2012) | 
A. I. Titov, H. Takabe, B. Kampfer , 
|Phys. Rev. Lett. 108, 240406 (2012)) 
T. Heinzl, A. Ilderton, and M. Marklund, 
|Phys. Lett. B 692, 250 (2010) | 

G. D. Tsakiris, K. Eidmann, J. M. ter Venn, and 

F. Krausz, |New J. Phys 8, 19 (2006)] 

G. Sansone, L. Poletto, and M. Nisoli, 
|Nature Photon. 5, 655 (2011)] 

J. H. Sung, S. K. Lee, T. J. Yu, T. M. Jeong, and J. Lee, 
|Opt. Lett. 35, 3021 (2010)1 

D. Herrmann, L. Veisz, R. Tautz, F. Tavella, K. Schmid, 
V. Pervak, and F. Krausz, |Opt. Lett. 34, 2459 (2009) | 

H. -C. W u, J. Meyer-ter Vehn, J. Fernandez, a nd B. M. 
Hegelich, |Phys. Rev. Lett. 104, 234801 (2010)] 

T. Z. Esirkepov, S. V. Bulanov, M. Kando, 

A. S. Pirozhkov, and A . G. Zhidkov, 

Phys. Rev. Lett. 103, 025002 (2009)| 



[14 
[15 
[16 

li- 
tis 

[19 
[20 
[21 

[22 
[2.3 

[24 

[25 
[26 
[27 
[28 
[29 
[30 
[31 
[32 



L. L. Ji, B. F. Shen, D. X. Li, D. Wang, Y. X. Leng, 
X. M. Zhang, M. Wen, W. P. Wang, J. C. Xu, and Y. H. 



Yu.. Phys. Rev. Lett. 105. 025001 (2010) 



L. L. Ji, B. F. Shen, X. M. Zhang, F. C. Wang, Z. Y. Jin, 

C. Q. Xia, M. Wen, W. P. Wang, J. C. Xu, and M. Y. 

Yu, |Phys. Rev. Lett. 103, 215005 (2009) | 

S. V. Bulanov, T. Esirkepov, and T. Tajima, 

|Phys. Rev. Lett. 91, 085001 (2003)| 

S. A. Skobelev, A. V. Kim, and O. Willi, 

|Phys. Rev. Lett. 108, 123904 (2012) [ 

A. Einstein, | Annalen der Physik 14, 194 (2005)] 

V. A. Vshivkov, N. M. Naumova, F. Pegoraro, and S. V. 

Bulanov, |Phys. Plasmas 5, 2727 (1998)| 

A. Macchi, S. Vcghini, T. V. Liseykina, and F. Pegoraro, 

|New J. Phys. 12, 045013 (2010)| 

P. Mulser and D. Bauer, High Power Laser-Matter In- 
teraction, Springer Tracts in Modern Physics, Vol. 238 
(Springer, 2010). 

T. Esirkepov, M. Borghesi, S. V. Bulanov, G. Mourou, 
and T. Tajima, |Phys. Rev. Lett. 92, 175003 (2004)] 
C. Thaury, F. Quere, J. -P. Geindre, A. Levy, T. Cec- 
cotti, P. Monot, M. Bougeard, F. Reau, P. d'Oliveira, 
P. Audebert, R. Marjoribanks, and P. Martin, 
|Nature Phys 3, 424 (2007)] 

C. Rodel, M. Heyer, M. Behmke, M. Kiibel, O. Jfickel, 
W. Ziegler, D. Ehrt, M. Kaluza, and G. Paulus, 
|Appl. Phys. B: Lasers Opt. 103 295 (2011) 



F. Pegoraro and S. V. Bulanov, 
|Phys. Rev. Lett 99, 065002 (2007)| 

AT V Krashcninnikov and F. Banhart, 
|Nature Mater . 6, 723 (2007)| 

W. van Dorp, C. Hagen, P. Grozier, B. van Someren, and 
P. Kruit, |Microelectron. Eng. 83, 1468 (2006)| 
M. Chen, A. Pukhov, T. P. Yu, and Z. M. Sheng, 
|Phys. Rev. Lett. 103, 024801 (2009")] 
M. Tamburini, F. Pegoraro, A. Di Piazza, C. H. Keitel, 
and A. Macchi, |New J. Phys. 12, 123005 (2010) | 
M. Tamburini, T. V. Liseykina, F. Pegoraro, 
A. Macchi, Phys. Rev. E 85, 016407 (2012)| 



and 



L. D. Landau and E. M. Lifshitz, "The classical theory 
of fields," (Elsevier, Oxford, 1975) Chap. 76, 2nd ed. 
M. Tamburini, F. Pegoraro, A. Di Piazza, 
C. H. Keitel, T. Liseykina, and A. Macchi, 



Nucl. Inst. Meth. Phys. Res. A 653, 181 (2011) 



